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We present two Fermi-Bose models with an approximate supersymmetry and 
which can be solved numerically with great accuracy using the renormalization 

C""- group method. The bosonic parts of these models consist in Dyson's hierar- 

chical model with one and two scalar components respectively. We discuss the 
question of the perturbative cancellations of divergences and compare with 
the non-perturbative fine-tunings necessary to keep the renormalized scalar 

(T) , mass small in cut-off units. We show evidence for non-perturbative cancel- 

lations of quantum corrections, however, we were not able to achieve exact 

(^ ' cancellations without fine-tuning. 
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The fact that the bare parameters of a scalar field theory require a fine-tuning in order to 
keep the renormalized mass small in cut-off units is usually regarded as an argument against 
fundamental scalars M. A possible resolution of this inelegant feature consists in adding 
degrees of freedom in such a way that the quantum fluctuations cancel, making small scalar 
masses a more natural outcome. 

From a non-perturbative renormalization group analysis, a relevant direction is nec- 
essary in order to describe massive particles, making the fine-tuning process unavoidable. 
This situation can be seen very simply in the case of a free scalar theory with a cut-off. In 
this example, the fact that the bare mass (mg) has to be small (in cut-off units) in order 
to get a small physical mass (ttir) is obvious since these quantities are identical. On the 
other hand in an interacting scalar theory, we usually need to take m 2 B negative and large 
in absolute value and also to adjust many digits of this quantity in order to get a small tur. 
The difference between the two situations is that in the second case, there is no bare quantity 
which controls the size of ttir. One would like to understand under which circumstances the 
inclusion of fermions allows us to obtain a small thr whenever we choose a small m B . 

There are known four- dimensional examples || where one can cancel the perturbative 
quadratic divergences by imposing simple relations between the Yukawa couplings and the 
scalar quartic couplings. However, it is not clear that there exists a non-perturbative regu- 
larization which fully preserves the perturbative naturalness. 

In the following, we present two models where there is an equal number of fermions and 
bosons and which can be solved non-perturbatively with great accuracy ||. The bosonic 
part of these model is a Dyson's hierarchical model J|]. In this model, the renormalization 
group transformation maps the local measure into another local measure. The price to pay 
for this simplifying feature is that the kinetic term is not ultra-local. The free action for TV 
massless scalar fields <{$ reads 

Si- = \ E <^Dly€ , (1) 

x,y,i 

where x and y run over the sites and % from 1 to N, The explicit form of D 2 xy is given below 
in Eq. (^|). The action for free massless fermions reads 

sfr = £ #D^« , (2) 

x,y,i 

where the ip® and $W are Grassmann numbers integrated with a measure 

n#F#i i} - (3) 

x,i 

As indicated by the notation, we have 

D 2 xy = Y,D xz D zy . (4) 

z 

The free action S B ee + S F ree is invariant at first order under the transformation 



S^=eJ2D xy ^ (5) 

a; 

02 = *Y, D **i>f' 



The e and e are Grassmann numbers. Integration by part or Leibnitz's rule cannot be used 
for D xy and the order ee variations do not cancel. 

We now give the explicit form of D 2 at finite volume. For a hierarchical model with 
2 Umax sites, we label the sites with n max indices x nmax , ...,X\, each index being or 1. In order 
to understand this notation, one can divide the 2 nmax sites into two blocks, each containing 
2 n max-i s it es# jf x nmax = 0, the site is in the first box, if x nmax = 1, the site is in the second 
box. Repeating this procedure n times (for the two boxes, their respective two sub-boxes, 
etc... ), we obtain an unambiguous labeling for each of the sites. With these notations, 

s'r = -p. Y. (f )" Y. ( Y. #g -o)" + jr| £ (4L.,,->) 2 

n=\ x nmax ,...,x n+1 ,i x n ,...,xi - -« x„ maa .,...,a:„ + i,j 

(6) 

The index n corresponds to the interaction of the total field in blocks of size 2 n . The constant 
cb = 2 1 ~ 2 ' D is a free parameter which controls the decay of the iterations with the size of 
the boxes and can be adjusted in order to mimic a D-dimensional model. Similarly the free 
massless fermionic action reads 

77"max ^ 

#" = -/fc£(f)" E ( E €_.-))( E 4L^)) (7) 

n=l E nmoa! ,...,a:„ + i,t £n,-,£i »«,— ,*i 

with cf = 2 1 ~ 1 / D . Using the techniques explained in ||, one can show that the fermionic 
operator is the square root of the bosonic operator (see Eq. @) provided that 



2 — cp 2 — c B 

We now introduce local interactions. The Grassmann nature of the fermionic fields 
restricts severely the type of interactions allowed. For instance, for one flavor (N = 1), the 
most general local measure is 

W(0,V^) = W(0) + v<M(0) (io) 

For convenience, we will always reabsorb the second term of Eqs. (||) and (^j) which are local, 
in the local measure. In the following calculations, W{<p) will take the Landau- Ginzburg 
(LG) form: 

W{<j>) oc e- (( ^f )+» 2+ ^ 4 ) . (11) 



If the two functions W and A are proportional, the fermionic degrees of freedom decouple. 
The renormalization group transformation takes the form 

W^2A*W (12) 

A^2[3 F A*W+(— )W*W (13) 

cf 

where the * operation means a convolution, a multiplication by an exponential and a rescal- 
ing of the new field. More precisely 

A * 5(0) = e^ 2 > / W 02Cb ^V( (02CB 3 + 0/) ) , (14) 

The introduction of a Yukawa coupling can be achieved by allowing a linear term in A(<f>). 
Such a term breaks explicitly the Z 2 symmetry of the LG measure. Such a model is charac- 
terized by a sudden change from the symmetric phase behavior to the broken phase behavior 
followed by unexpectedly long low-temperature "shoulders" (see M for an explanation of 
this terminology). 

A richer behavior is observed in the case of the two flavors (i = 1,2) models. In the 
following we have restricted our investigation to the type of bilinear coupling appearing in 
the Wess-Zumino [|J model, namely 

w(<£ (i) ,v (0 ,$ (i) ) = w(0 (i) ) + a(0 w )(^ (1 V (1) + v> (2 V (2) ) + 

5((/)«)^ (1 V (2) - b*(^ w ) W (2) + t(0«)^ (1 V (1) V' (2 V (2) • (is) 

For convenience, we again absorb the local parts of Eqs. (||) and (§). This is not the most 
general measure, however it closes under the renormalization group transformation which 
takes the form 

W -> (W *T + A* A + B * B*) = W' 

A -> p F W' + —A * T 
cf 

4 , , 

B^—B*T (16) 

cf 

T ^^T*T + (3f—A *T+ {(3 f ) 2 W . 
cp c F 

In addition we will impose that the function B have the following form: 

B{^) = (0« + z</>( 2 ))P((0 (1) ) 2 + (0 (2) ) 2 ) , (17) 

while W, A and T are 0(2)-invariant. The model is then invariant under the R-symmetry 

(0(l) H _ i 0(2))_^ e ^ (0 (l) H _ # (2) ) 

^U) _+ e -*fytf) (18) 



We now present three numerical calculations performed with the second model. In all 
cases we will set D = 4 in c B and cp. In the following, we have chosen the value of (3b and 
(3p in such a way that 

(3 F c F f f3 B c B \i _ -. /-|q\ 



l — Cp £ — c B 

in order to make the perturbative expansion more similar to usual Feynman diagram's 
calculations. 

First, we consider the case where the fermions decouple from the bosons. W takes a LG 
form 

W(CJ>) OC e -«^) + ^l)E i (^) 2 + A B (E,(* W ) 2 ) 2 ) . (2Q) 

The value of m R , defined as the inverse of the zero-momentum two-point function, is shown 
in Fig. 1. as a function of m B . These quantities are expressed in cut-off units. For reference 
we have also displayed the one-loop perturbative result and the trivial gaussian result. One 
sees that the scalar self-interaction moves m R up and m R ~ 0.2 when m B goes to zero. The 
one-loop result is quite good when m R is large enough but deteriorates when this quantity 
becomes smaller. 

In the second calculation, we consider a bosonic model with a bare mass mg and Xb = 
coupled to a fermion with the following couplings: 

A= (-l-m B )W 

P = g y W (21) 

T = ((-1 - m B f + ^((0 (1) ) 2 + (<P (2) f))W . 



The results are shown in Fig. 2 for g y = V0.08 ~ 0.28. One sees that the Yukawa coupling 
moves m R down. For m B — 0.094, m R becomes and for smaller of m B , we enter the broken 
symmetry phase. 

We have then repeated the second calculation with X B = 0.01 instead of 0. In perturba- 
tion theory, the one-loop quadratic divergence cancel when m# = and 

8Xb = g 2 y , (22) 

which justifies our choice of coupling constant. The results are shown in Fig. 3. One sees 
that the Yukawa coupling in part cancels the effects of the scalar self-interaction, however, 
the cancellation is not as good as in the one-loop formula where m R goes to zero when m B 
goes to zero. Instead, we found numerically that m R ~ 0.044 when m B goes to zero. A 
summary of the three numerical results is shown in Fig. 4. 

It is possible to fine-tune g y in order to get m R = 0. An example is shown in Fig. 
5 for Xb = 0.01 and m 2 B = 0.01. We see that there exist a critical value of g y which is 
approximately 0.46 and where m R becomes 0. For larger values of g y , we enter the symmetry 
broken phase. An essentially similar figure is obtained for m 2 B = 0. In both cases, the exact 
critical value of g y is about 50 percent larger than the perturbative one. 

In conclusion, we have shown that the idea of canceling the quantum correction inspired 
by perturbation theory have qualitatively a non-perturbative counterpart. However, we have 
not found a way to make this cancellation very accurate or exact without fine-tuning. 
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FIG. 1. The renormalized mass as a function of the bare mass in a bosonic 0(2) model with 
bare quartic coupling fixed to 0.01. 
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FIG. 2. The renormalized mass as a function of the bare mass in a bosonic 0(2) model with 
bare quartic coupling fixed to and a Yukawa coupling equal to \/0.08 
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FIG. 3. The renormalized mass as a function of the bare mass in a bosonic 0(2) model with 
bare quartic coupling fixed to 0.01 and a Yukawa coupling equal to \/0.08. 
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FIG. 4. The numerical results of Figs. 1, 2 and 3 combined together. 
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FIG. 5. The renormalized mass as a function of the Yukawa coupling in a bosonic 0(2) model 
with bare quartic coupling fixed to 0.01 and a bare mass m 2 B = 0.01. 
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